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ABSTRACT 



By using the geostrophic approximation in various of the horizontal 
acceleration terms in the Eulerian equations for horizontal flow, wind 
components are computed from which relative vorticity is determined. 

The problem of obtaining a stream function is then reduced to solving 
a Poisson-type equation, thus avoiding the use of the balance equation. 

Two different computation schemes for the x, y components of the 
500-mb wind are investigated. The stream functions obtained from these 
methods are then compared with those obtained through the solution of 
the balance equation. The results of forecasts for 24 and 48 hours 
using the various stream functions with the barotropic model are also pre- 
sented. 

The writer wishes to express her appreciation to Professor George 
J. Haltiner of the U. S. Naval Postgraduate School for his assistance 
and guidance in this investigation. 

Appreciation is also expressed to the personnel of the U. S. Navy 
Fleet Numerical Weather Facility for their cooperation during the prepa- 
ration of this paper. The writer is especially indebted to Mr. Geirmun- 
dur 'Arnason for his many informative discussions and suggestions, and 
to Mr. Milton H. Reese for his assistance in programming. 
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the Coriolis parameter,, 3- -O- Sun 0, where <$> is the 
geographical latitude. 

the mean value of the Coriolis parameter, 3l jx. sun 45°* 

the height of an isobaric surface. 

the upward component of the apparent gravitational 

acceleration. 

the geopotential. • 

the relative vorticity. §= ^ < 

the absolute vorticity. ^+-f ♦ 

the stream function for the non-divergent component of 
velocity. 

horizontal Laplacian operator on a constant pressure 
surface . 

horizontal Jacobian operator, -J C\B) = ^ ~ 
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1. Introduction. 



It has been found that the use of the geostrophic approximation in 
the quasi-geostrophic barotropic model for numerical prediction gives 



such errors can be eliminated by replacing geostrophic wind and geostro- 
phic vorticity by wind and vorticity derived from a stream function, thus 
retaining only the divergence-free part of the wind 

The present method used by both the Joint Numerical Weather Predic- 
tion Unit (JNWPU) in Suitland, Md. , and the U. S. Navy Fleet Numerical 
Weather Facility (FNWF) , Monterey, Calif. , to obtain a stream function 
is through the solution of the balance equation 



where $ is the geopotential. This equation, a non-linear partial differen- 
tial equation of the Monge-Ampere type, may be elliptic in one part of the 
map, parabolic in another and hyperbolic in a third area. The problems 



encountered in obtaining significant solutions to this equation have been 



The purpose of this research is to examine the possibility of obtaining 
a stream function in a somewhat simpler manner. By using the geostro- 
phic approximation in various of the horizontal acceleration terms in the 
Eulerian equations for horizontal flow, wind components may be computed 
from which the relative vorticity may be determined. The problem of ob- 
taining the stream function is then reduced to solving the equation 



rise to spurious anticyclogenesis [l] . It has further been shown that 




discussed by 'Arnason 







( 2 ) 



a Poisson-type equation readily solved by numerical methods. Thus the 
use of the balance equation may be avoided. 
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2. Background. 



Assuming the hydrostatic approximation, and using pressure (p) as 
the vertical coordinate, the equation of motion may be expressed as 



In scalar component form, equation (3) becomes, upon expanding 
the acceleration components , 



Various approximations have been used to represent the actual wind 
The simplest is the geostrophic approximation in which the entire accele 
ration is assumed to be identically zero, yielding 



A second approximation widely used is that the flow is gradient. 
With this assumption, tangential accelerations are omitted, but centri- 
petal accelerations are retained. 

The wind law to be used herein primarily consists of omitting the 
local time derivatives. In addition, since the results are to be applied 
to one level only in the present study, the acceleration terms involving 
the vertical velocity advection will be omitted. Thus equations (4) 



d^=-3%z-f IkxV . 



(3) 




(4) 



where uJ - • 






(5) 



reduce to 








( 6 ) 



ax ^ 
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By introducing geostrophic wind components in various terms on the 
left-hand side of equation (6) the u, v components of the wind and the 
relative vorticity were obtained. 

Once the relative vorticity had been computed , a stream function, 
y , was obtained by solving the Poisson-type equation 

S = S- (7) 

g 

using the relaxation technique. The introduction of the coefficient -r~ 

f 

is necessary to give the computed y-field units of height. This vj;-field 
as computed was then used to make 24- and 48-hour barotropic prognoses 
which were compared to forecasts of both the JNWPU and FNWF stream 
functions using the same barotropic model. 

All computations were done on the Control Data Corporation Model 
1604 Computer using the U. S. Navy Fleet Numerical Weather Facility's 
octagonal grid. This grid, centered on the pole and extending to approxi- 
mately 10° N latitude, consists of 1977 grid points with a grid size of 
381 km at 60° latitude. Geostrophic winds and geostrophic vorticity 
were also computed in order that the effects of the approximations intro- 
duced could be compared. Since it is known that the geostrophic approxi- 
mation is poor in low latitudes, and that small errors in the pressure field 
will give relatively large errors in computed winds in this area, the sine 
of the latitude was not allowed to become smaller than sine 15° through- 
out the computations . 

Initial data used were the 500-mb heights for OOZ„ 3 January 1958. 
The v^-fields for OGZ, 4 and 5 January 1958 were also computed in order 
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to compare the forecast fields with the computed. 
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3. Procedures. 



Basically, two computation schemes were investigated to obtain 
approximations of the x, y components of the 500-mb wind. In the first 
method, geos trophic wind components were substituted in the partial de- 



rivatives of equations (6) . Thus the following simultaneous system of 
equations is obtained: 



These equations were then solved for u and v, with the following re- 
strictions imposed: 

(a) an arbitrary lower limit of .25 was imposed on the terms 



inertial instability; 

and (b) the determinant of the coefficients was not allowed to become 
smaller than .125, to avoid excessive wind speeds. The above method 
will hereinafter be designated as Method I. 

The second method investigated, to be referred to as Method II, 
consisted of replacing all the velocity components in the horizontal ac- 
celeration terms of equations (6) by the geostrophic components. Equa- 
tions (6) then become 



This method, though presumably a cruder approximation to the actual 




( 8 ) 





(9) 
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wind, has the decided advantage of avoiding the arbitrary restrictions 
imposed in Method I, as well as a greater simplicity in programming. 

Since the stream field obtained by Method II was very nearly the 
same as that obtained by Method I , one can consider the sequence of 
approximations obtained by the equations 

(n) 

U =Uo- £ *2- - f ^ 

3 f f 

(10) 



(r») (in 

V ” Vj +■ 



\ (n-0 

Jv 


V dv 


^x 


f ^ 


V (*•>) 


V dU 


^X 


f 



(x) (a) 

Method III consisted of computing u and v 



with h = t * • - • - 

After computing the wind fields for the various methods, and the 
geostrophic wind,\y^ - - Jjf ^,2 XIK, for OOZ, 3 January 1958,- the 
difference between computed and geostrophic wind speeds, I Vl ~ l^gl, 
was examined in order to assess the effects of the various approximations 
on the horizontal acceleration. In the regions of sharply defined cyclonic 
curvature, Method II produced the greatest decrease, and, in regions of 
anticyclonic curvature, the least increase in the geostrophic wind. This 
is to be expected as may be seen from gradient-wind considerations. 
When the geostrophic approximation is used in the acceleration terms, 
the latter will tend to be overestimated in regions of cyclonic curvature 
and thus reduce the wind too much below the geostrophic. Similarly, in 
the vicinity of a ridge, the geostrophic approximation in the acceleration 
will tend to underestimate the acceleration and thus the resulting in- 
crease of the wind over the geostrophic will be inadequate. 

The introduction of the higher order approximation by Method III 
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tended to bring the wind speed again toward the geostrophic and, in 
regions of extreme cyclonic curvature, the computed speeds exceeded 
the geostrophic in some cases. This will be further discussed in con- 
nection with vorticity. 

Tables 1 and 2 show representative values of the difference be- 
tween computed and geostrophic wind speeds found in regions of strong 
curvature . 
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Table I: Maximum difference between computed and geostrophic wind 

speeds observed in selected areas of well-defined cyclonic 
curvature 







Geostrophic 


Maximum reduction in 
geostrophic wind speed 




Location 




wind speed 


Method I 


Method II 


Method III 


Eastern 
United States 


50 to 58 mps 


-17 mps 


-26 mps 


-8 mps 


Eastern 

Atlantic 




30 to 38 mps 


-12 mps 


-19 mps 


-5 mps 


Eastern 

Europe 




30 to 38 mps 


- 8 mps 


-11 mps 


-4 mps 


Off East 
Coast of Asia 


50 to 62 mps 


-20 mps 


-31 mps 


-9 mps 


Central 

Pacific 




50 to 60 mps 


-19 mps 


-28 mps 


-7 mps 


Table II: 


Maximum difference between computed and geostrophic wind 
speeds observed in selected areas of well-defined anti- 
cyclonic curvature 






Geostrophic 


Maximum increase in 
geostrophic wind speed 




Location 




wind speed 


Method I 


Method II 


Method III 


Off East Coast 
of North America 


15 to 20 mps 


13 mps 


7 mps 


9 mps 


Siberia 




25 to 31 mps 


8 mps 


6 mps 


7 mps 


Central 

Pacific 




30 to 39 mps 


18 mps 


11 mps 


12 mps 


Western 

Canada 




4 to 14 mps 


3 mps 


2 mps 


2 mps 
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Next the fields of geostrophic vorticity, —• <7 Z , and relative 

^ T 

vorticity,, $= , were computed. As was expected, the de- 

rived values of cyclonic vorticity associated with the migratory lows and 
troughs were smaller than the corresponding geostrophic quantities, with 
the reverse true in the case of anticyclonic vorticity. As shown in Table 
III, Method II produced the greatest decrease in cyclonic vorticity and 
the least increase in anticyclonic vorticity. With the introduction of the 
higher-order approximation, Method III, the values of both the positive 
and negative vorticity centers increased over those obtained from the 
first two methods. Again, in regions of strong cyclonic shear, the maxi- 
mum values exceeded the geostrophic. Further iterations using the tech- 
nique of Method III were made to determine if the introduction of still 
higher-order approximations would indicate convergence. 

The vorticity fields computed from u <J and v ' 5) did show some reduc- 
tion in the cyclonic vorticity maxima, but the negative vorticity values 
kept increasing. With the next iteration, u^' and v (,f) , the values ob- 
tained in regions of strong relative vorticity produced overflow in the 

fixed-point computations of this parameter, which had been scaled for a 

-6 -1 

maximum value of 253 x 10 sec . This suggests, perhaps, that the 
sequence of computed winds was not converging toward realistic values, 
at least in one particular region, having given rise to excessive vorti- 
city values . 

After the relative vorticity fields were computed, the Poisson equa- 

. c U. 

tion j = V W-* was solved for the stream function using the 
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relaxation technique. This technique as programmed by FNWF scans 
through the field of (JJ laterally using the extrapolated Liebmann method 
of computing a new point value 









(H) 



where Rij is equal to ^ (>Ap -If.i) and is the relaxation factor, 

4 j J 4 

in this case 0 » 35 . The 500-mb height field was used as the initial ap- 



proximation. 
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Table III 



Position of 
associated low 


s 

°3 


Maximum value of center of 
cyclonic relative vorticity 
Method I Method II 


Method III 


39N 


10E 


67 


52 


47 


60 


58N 


22E 


80 


68 


63 


70 


38. 5N 


33E 


80 


63 


58 


73 


78. 5N 


88E 


38 


34 


34 


36 


54N 


147E 


73 


61 


57 


65 


56N 


164W 


59 


51 


50 


51 


84N 


135W 


51 


51 


50 


52 


55N 


71W 


84 


75 


71 


74 


46N 


23W 


100 


71 


59 


92 
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4. Results of Stream-Function Computations. 



A comparison of the computed stream-function fields with those ob- 
tained by the FNWF solution of the balance equation reveals that, in all 
instances, the gradient of y between the pole and the grid boundaries is 
slightly greater for the fields computed by the methods of this study than 
for those computed through the balance equation, the difference varying 
from 200 to 400 feet in the cases studied Through the relationship 
xiK , it is obvious, then, that on the average the non- 
divergent part of the wind, computed by the methods described herein, 
slightly exceeds that obtained through the balance equation. Since, 
however, the difference in the gradient is almost uniform from the bound- 
aries to the pole, the difference in wind speed at any point would be only 
of the order of centimeters per second. 

Considering only the stream-function fields computed by the methods 
of this study, the gradient of vy between pole and grid boundaries was 
least for Method II, corresponding to the lowest wind speeds. 

A comparison of the computed stream functions with the FNWF 
stream functions and the height fields themselves reveals that the place- 
ment of the major systems is the same in all cases. However, an exam- 
ination of the stream fields obtained by Method I reveals that the decrease 
in the amplitude of the troughs is slightly greater than that observed be- 
tween the height and the FNWF stream-function fields. This decrease 
in amplitude of the major troughs is extremely pronounced in the 
vy -fields computed by Method II, presumably due to the underestimation 
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of the wind by this method as previously mentioned. 



With the introduction of the higher order approximation, Method III, 
this pronounced weakening of the troughs is not observed, the amplitude 
of the features being approximately the same as those of the FNWF 
vp -fields . 
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5. Results of Stream Barotropic Prognoses. 



Using the stream functions obtained by the various methods, 24- 
and 48 -hour prognoses were made using the barotropic forecasting equa- 
tion 






( 12 ) 



where rJ~ is the Jacobian operator [3j . The term - is added to 

the non-divergent barotropic model to stabilize the ultra-long atmos- 
pheric waves. In these computations x was assigned the value of 4, 
in accordance with the results of Cressman L4] • 

In order to have some basis for comparison, 24- and 48-hour prog- 
noses were also made using the same model with the stream functions 
obtained through both the FNWF and the JNWPU solutions of the balance 
equation. For verification purposes, difference fields were formed and 
the "pillow" and root-mean-square error computed over the entire grid 
using the formulae 



Pillow = 



isit, 

Z (a-B)„ 






RMSE 




(13) 

(14) 



where A is the computed stream function value and B is the forecast 
value . 

As can be seen from Table IV, no significant difference in either 
the "pillow" or root-mean-square error is apparent for the case studied. 
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However, it may be noted that the root-mean-square error for both the 



24- and 48 -hour forecasts is least for Method II and largest for Method 

III. 

An examination of the error fields themselves showed a marked simi- 
larity in the regions of errors greater than 200 feet» The only major dif- 
ferences that could be noted were in an area north of the Hawaiian 
Islands, where Method III had a much larger error than other methods, 
and in an area south of Greenland, where both Methods I and II produced 
errors approximately 200 feet less than the other methods. 

Since the barotropic model used for forecasting does not include 
provisions for development, only the movement of the major systems was 
used to compare the various methods. Considering only this aspect, no 
great differences appeared with regard to most of the major systems . The 
only notable exceptions were the troughs in the Eastern Pacific and the 
Eastern Atlantic where movement was overestimated less using the 
FNWF stream function than by the other methods. In addition, there was 
a slight indication that, at least in this case, the systems tended to 
move slightly faster when Method I was used. 
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Table IV 



Method used to obtain 
initial stream function 


24-hour baro- 
tropic forecast 


48 -hour baro- 
tropic forecast 




Pillow 


RMSE 


Pillow 


RMSE 


FNWF solution of 










balance equation 


+ 8 ft. 


174 ft. 


+ 6 ft. 


272 ft. 


JNWPU solution of 










balance equation 


-36 ft. 


172 ft. 


-70 ft. 


273 ft. 


Method I 


-14 ft. 


169 ft. 


-56 ft. 


282 ft. 


Method II 


+9 ft. 


164 ft. 


+24 ft. 


254 ft. 


Method III 


-17 ft. 


182 ft. 


-62 ft. 


311 ft. 
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6. Conclusions. 



Since the methods described herein have been applied to only one 
synoptic situation in this study, only tentative conclusions can be made 
at this point as to the relative merits of the various ways of obtaining a 
stream function for use in numerical prognostic models. All methods 
have shown comparable results and, until many different weather pat- 
terns are tested,, it is difficult to separate the errors in the forecasts 
into those arising from inadequacies in the barotropic model and those 
due to the assumptions made in the various methods. 

It has been seen, however, that the stream-function fields obtain- 
ed from Method II had the least computed root-mean-square error of any 
of the methods tried and a "pillow" smaller than those computed using 
the other two methods of this study or the JNWPU solution to the balance 
equation for both the 24- and the 48-hour prognoses. However, a greater 
decrease in amplitude in the transition from the height to the stream- 
function field has also been observed when this method was used and 
must also be considered in evaluating this method. Although this method 
presumably underestimates the total wind speed, in the case studied the 
24-hour prognostic movement of the major systems did not differ appre- 
ciably from that of the other methods used. Since no arbitrary restric- 
tions need be imposed when this method is used, and since it is easier 
to program , further study is certainly warranted. 

Since the introduction of the higher-order approximation into the 
same basic computation scheme (Method III) produced spurious results 
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in both the wind and the relative vorticity fields, the use of this pro- 



cedure appears questionable . However, the introduction of certain re- 



in Method I, may possibly eliminate the difficulties encountered 

The procedure described under Method I apparently has some merits, 
although further study should be made of the arbitrary restrictions placed 
on this method for the computations of this study. It would also be most 
interesting to try further iterations of this method using equations of the 



to determine what the effect on the geostrophic deviation and relative 
vorticity would be in comparison to the effects observed upon introduc- 
tion of the higher approximations in Method III. Since extensive repro- 
gramming would be involved in doing this, such iterations were not 
attempted in this study. 



\ («) 

strictions on the values of 4^- and 




form 





(15) 
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Appendix I 



The following charts were computed during the course of this in- 
vestigation and are available upon request from the U. S. Naval Post- 
graduate School, Monterey, California. 

(1) 500-mb heights, OOZ, 3 Jan. 1958 

(2) 500-mb heights , OOZ, 4 Jan. 1958 

(3) 500-mb heights , OOZ, 5 Jan. 1958 

(4) Geostrophic Wind Speed, OOZ, 3 Jan. 1958 

(5) Method I, 500-mb Wind Speed, OOZ, 3 Jan. 1958 

(6) Method II, 500-mb Wind Speed, OOZ, 3 Jan. 1958 

i 

(7) Method III, 500-mb Wind Speed, OOZ, 3 Jan. 1958 

(8) Geostrophic Relative Vorticity, OOZ, 3 Jan. 1958 

(9) Method I, Relative Vorticity, OOZ, 3 Jan. 1958 

(10) Method II, Relative Vorticity, OOZ, 3 Jan. 1958 

(11) Method III, Relative Vorticity, OOZ, 3 Jan. 1958 

(12) FNWF Stream Function, OOZ, 3 Jan. 1958 

(13) Method I, Stream Function, OOZ 3 Jan. 1958 

(14) Method II, Stream Function, OOZ, 3 Jan. 1958 

(15) Method III, Stream Function, OOZ, 3 Jan. 1958 

(16) FNWF Stream Function, OOZ, 4 Jan. 1958 

(17) Method I, Stream Function, OOZ, 4 Jan. 1958 

(18) Method II, Stream Function, OOZ, 4 Jan. 1958 

(19) Method III, Stream Function, OOZ, 4 Jan. 1958 

(20) FNWF Stream Function, OOZ, 5 Jan. 1958 
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(21) Method I, Stream Function, OOZ, 5 Jan. 1958 

(22) Method II, Stream Function, OOZ, 5 Jan. 1958 

(23) Method III, Stream Function, OOZ, 5 Jan. 1958 

(24) FNWF 24-hr. Stream-Function Prognosis, Verifying Time, 
OOZ, 4 Jan. 1958 

(25) Method I, 24-hr. Stream-Function Prognosis , Verifying 
Time, OOZ, 4 Jan. 1958 

(26) Method II, 24-hr. Stream-Function Prognosis, Verifying 
Time, OOZ, 4 Jan 1958 

(27) Method III, 24-hr. Stream-Function Prognosis , Verifying 
Time, OOZ, 4 Jan. 1958 

(28) Difference between FNWF Computed ^ “field and 
Prognosis, OOZ, 4 Jan. 1958. 

(29) Difference between JNWPU Computed -field and 
Prognosis, OOZ, 4 Jan. 1958 

(30) Difference between Method I Computed -field and 
Prognosis, OOZ, 4 Jan. 1958 

(31) Difference between Method II Computed ^ -field and 
Prognosis, OOZ, 4 Jan. 1958 

(32) Difference between Method III Computed vy -field and 
Prognosis, OOZ, 4 Jan. 1958 
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